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TECHNICAL NOTE 5057 


A SIMPLIFIED MATHEMATICAL MODEL FOR 
CALCULATING AERODYNAMIC LOADING AND DOWNWASH 
FOR WING-FUSELAGE COMBINATIONS WITH WINGS 
OF ARBITRARY PLAN FORM^ 

By Martin Zlotnick and Samuel W. RoLinson, Jr. 

SUMMARY 


For the pui^jose of calculating the aerodynamic loading on the 
fuselage, the midwing wing-fuselage combination with a fuselage of cir- 
cular cross section can he represented hy a sin 5 )le system of horseshoe 
vortices located on the wing with images located inside the fuselage. 

By using this sin^lified mathematical model or the extension of it given 
in an appendix for nonmidwing configurations with fuselages of arhitrary 
cross section, a method for ceilculating the lift and longitudinal center 
of pressure on the fuselage in the presence of the wing at subsonic 
speeds is presented. 

In addition the report shows how the siniplified mathematical model 
can he used for calculating the downwash behind the wing and for calcu- 
lating the spanwise lift distribution on the wing for midwing configura- 
tions with axisymmetric fuselages. 


INTRODUCTION 


Mitual Interference between wing and fuselage has a significant 
effect on the pitching moment of the wing-fuselage combination, since 
the longitudinal distribution of the aerodynamic loading on the fuse- 
lage is altered by the presence of the wing. Multhopp (ref. l) has 
developed a theoretical method for calculating the pitching moment on 
wing-fuselage combi nations with unswept wings which gives good sigree- 
ment with experimental results. However, for calculating the pitching 
moment on swept-wing configurations, only semien^lrical methods such 
as that of reference 2 are available. 

^Supersedes the recently declassified NACA RM L52J27a, "a Sinpli- 
fled Mathematical Model for Calculating Aerodynamic Loading and Downwash 
for Midwing Wing-Fuselage Conibinations With Wings of Arbitrary Plan Form" 
by Martin Zlotnick and Samuel W. Robinson, Jr., 1955- 
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The aerodynamic loa d ing on the fuselage in subsonic flow is con- 
sidered in references 1 and 2 to he made up of two parts. The first 
part, which is due to the fuselage angle of attack resulting from both 
the geometric angle of attack of the fuselage and the upwash angle 
induced by the wing, has been calculated in reference 1 for unswept- 
wing configurations and modified in reference 2 for the swept-wing 
case. The second part, which is often referred to as the wing lift 
"carried over" by the fuselage (and which will be referred to in this 
paper as the "induced lift"), has been calculated analytically for the 
Tinswept-wing case in reference 1. However, for swept wings this calcu- 
lation cannot be applied and no other theoretical method has been avail- 
able. In reference 2 this component of loading is estimated by an 
en^iirical method. 

In order to calculate the induced lift on the fuselage in combina- 
tion with a wing of arbitrary plan form, a method is suggested in the 
present paper which is based on Lemiertz' theoretical work (ref. 3)» 
Lennertz^ results, which are concerned only with an unswept lifting 
line passing through the axis of an infinitely long cylindrical fuselage, 
are, in effect, generalized so that for midwing configurations in sub- 
sonic flow the magnitude, lateral distribution, and longitudinal center 
of pressure of the induced lift on an axisymmetric fuselage combined 
with a swept lifting line, or even a lifting sturface, can be calculated. 
An extension of this method to configurations with arbitrary fuselage 
cross section and wing location is described in appendix A. This exten- 
sion is based on a resiilt of Flex (ref. 1|-) and permits the calculation 
of the magnitude and longitudinal center of pressure of the induced 
lift, but not its lateral distribution. A nxjmeidcal example is given 
in appendix B to Illustrate the induced-lift calculations. The effect 
of finite fuselage length is estimated qualitatively from results of an 
approximate calculation of the variation with fineness ratio of the 
induced lift on an ellipsoid of revolution combined with an infinite 
vortex. The approximate calculation is shown to give res\ilts which 
agree with results of Vandrey (ref. 5)* 

In calcTilating the induced 3J.ft on the infinitely long fuselage 
with a circular cross section, the wing-fuselage combination is replaced 
with a single system of horseshoe vortices (or doublets) on the wing, 
with images inside the fuseleige. This representation may be used, with 
some modifications, in calculating the lift on the wing for a midwing 
configuration with a fuselage of circular cross section and also the 
corresponding downwash. However, for calculating the downwash and the 
lift on the wing, the sin^jlified representation of the wing-fuselage 
combination, although it is considered to be adequate, is no longer 
rigorous, and the general applicability will depend on e^erimental 
verification. 
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SYMBOIS 


A 

a 

lo 

c 

c 

cz 

d 

Lf 

m 

P 

Po 

a 

s 

s 

Uo 

^max 

Xcp 

r 

p 


aspect ratio (wing alone) 

maximum radius of body of revolution 

wing span 

wing chord 

mean chord, S/b 

local lift coefficient 

major axis of ellipsoid of revolution 

total, lift on fuselage 

number of horseshoe vortices 

static press\ire 

free-stream static pressure 

dynamic pressure 

wing area (wing alone) 

semispan of bound leg of horseshoe vortex 
free-stream velocity 

maximum longitudinal, velocity on surface of body 
longitTidlnal center of pressure 
vortex strength 
mans density of air 


x,y,z longitudinal, lateral, and vertical ordinates, respectively 

u,v,w longitudinal, lateral, and vertical velocities, respectively 

zjj hei^t of plane of wing above fuselage axis 

Subscripts: 

a downwash point 


I 


on lower surface 
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n particular pair of horseshoe vortices 

u on upper surface 

Superscripts : 

' dimensionless with respect to a 

* dimensionless with respect to h/2 

A har over a symbol denotes that the quantity is dimensionless 
with respect to s. 


BASIC CONSIDERATIOKS 

Midwing Configurations With Axlsymmetric Fuselages 


In the major part of the analysis the fuselage is assumed to he an 
infinite cylinder. To obtain a qual itative estimate of the effect of 
the finite fuselage length, the variation with fineness ratio of the 
lift on an ellipsoid of revolution with a wing of infinite span will he 
calculated approximately in a subsequent section. An approximate method 
for making the smaT 1 correction for this effect will he indicated. 

Review of results obtained by Lennertz .- Lennertz (ref. 5 ) has 
calculated the lateral and longitudinal lift distribution on the fuse- 
lage of an idealized wing-fuselage configuration, in which the fuselage 
is represented by an infinite circular cylinder and the wing by a vortex 
having constant spanwlse circulation. The vortices trailing from the 
wing tips have images inside the cylinder, and the bound vortex is 
extended Inside the cylinder to join the trailing image vortices as 
shown in figture 1. With this configuration, the boundary condition of 
zero velocity normal to the surface of the cylinder is satisfied only 
at infinity and in the plane normal to the cylinder axis which passes 
through the vortex, so that it is necessary to superpose an additional 
potential, which is calculated in reference 5- For this case the lat- 
eral lift distribution is obtained by considering the momentum change 
in a vertical plane infinitely far behind the wing, and the longitudinal 
lift distribution is obtained by the use of Bernoulli's equation. 

Induced lift on infinite cylinder with finite wing having constant 
spanwise distribution of circulation .- In this section the fuselage 
lift I^ and its lateral distribution dLf/dy are shown to be xinaffected 
by the additional potential, so that only the con 5 >onents of Lf and 
dlxf/dy due to the vortex potential need be calculated. 
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With all other singularities neglected, the lift on a longitudinal 
section of the cylinder dy (see fig- l) due to the vortex potential 
is calculated as follows: 

The pressure p of any point on the surface of the cylinder can 
he written 


P = Po + SP 


¥ 


I^Uq + Au)^ + Av^ + Aw^ 


where Au, Av, and Aw are, respectively, the longitudinal, lateral, 
and vertical con^ionents of velocity induced by the vortices. The sec- 
tion lift dLf/dy is then written 


dl^ 


poo 

/ (Pu - 


where p^^ and p^ are the pressures on the upper and lower surfaces, 
respectively. Then 



Au dx 


( 1 ) 


Only the velocities induced by the bound vortex contribute to the lift. 


It may be noted that, since the distribution of the longitudinal 
velocities induced by the bound vortex is symmetrical about its axis 
at every section dy, the longitudinal center of pressure of the lift 
due to the vortex potential is on the axis of the bound vortex. 


A closed expression for the integral in equation (l) may be readily 
derived as follows : Consider the rectangular path indicated by the 

dashed line in the upper ri^t sketch in figure 1 and the cross section 
downstream at infinity shown in the lower right sketch of the same fig- 


ure. 


The line integral 




Au ds 


of the tangential velocity con^jonent 


teiken around the complete rectangle must be P where the path links 


one of the horseshoe vortices 


/|a2 


VI 


b/2 


< |y| < la|), and must be zero 


where it does not link one of the horseshoe vortices 0 < |y| < 


b/2 )■ 


Then the desired integral in equation (l), which is the longitudinal 
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portion of the cou 5 )lete line integral 


ds. 


most equal ^Au ds 


f 


minus the line integral along the short vertical line at infinity 
I 

Aw(t»,y,z)dz, vhere u refers to the upper surface and I refers 
u 

to the lower siirface. That is. 


X 


/“ n pZ 

An dx = (^Au ds - / Aw(“,y,z)dz 


(2) 


Since 


then 


oz 


f dz = - 01 

^u 


where 0^ antS 0^ gire the potentials on the upper and lower stirfaces, 


respectively, in the plane x = ». When 


h/2 


< |y| < |aU 


0u = -0Z 


= i4an-l - tan-1 'feg - 7^, . tan'l + tan'l 

a^ 


_ k 


y + I 


y - 


h/2 


y + 


= lrtan-l^i , 

/h^2 


\/a2 - y2 


2 


and 




Au ds = r 
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and when 0 < | y | < 
a2 


h/2 


h/2 

is increased by it: 


the strength of the Image vortex located at 


01 


X 

2it 


-i 1^2 _ y2 -1 \/a^2 _ y2 / , 1/3,2 _ y2\ -i \/a2 — 

tan'l y - - tan-1 ^ - (it + tan“l + tan"l 1 ^ 


y-| 


y.| 


y- 


h/2 


y+ 


h/2 


i r 

It 


ft an" 


-1 


2| /a2 - y2 


( 


If -a. 


'2 


and 




ds = 0 


Then, from eq.uatlon (2), the lateral lift distribution of the induced 
lift can be written 


dLf 


1 - — tan' 

Jt 


. 1 ^/ 


a2 - y2 


ir - 


(-a < y < a) ( 3 ) 


and after integration over y. 


Lf = 


2pUoT a 



(h) 


These expressions are the same, except for the difference in notation, 
as the expressions obtained in reference 3 for dl^/dy and the total 
fuselage lift Lf. The con^ionents of Lf and dl^/dy due to the addi- 
tional potential must therefore be zero. 
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The longitudinal lift distrihution on the cylinder calculated in 
reference 5^ which includes the effect of the additional potential as 
well as the vortex potential, is sjLlghtly different from that which 
would he calculated by taking into account the effect of the vortex 
potential alone. However, both distributions are symmetriceil about 
the axis of the bound vortex. The longitudinal center of pressure of 
the lift due to the vortex potential must be at the bound-vortex axis 
becatise the longitudinal velocities induced on the surface of the 
cylinder by the bound vortex are the same in front and in back of the 
bound vortex, as noted previously. That the longitudinal center of 
pressure of the lift due to the vortex potential and the additional 
potential must be at the axis of the bound vortex is indicated by the 
calculations of reference 3 (althou^ not explicitly stated) and can be 
shown as follows: 

Superposition of inf ini te vortices canceling the semi-infinite 
trailing vortices of figure 2(a) as shown in figure 2(b) will not change 
the longitudinal lift distribution on the cylinder, since no longitudinal 
velocities are induced. It is apparent that the system is the same as 
before, with the direction of the trailing legs of the horseshoe vortices 
reversed; therefore, the longitudinal distribution on the infinite 
cylinder must be symmetrical about the axis of the bound vortex, since, 
if it were not, the longitudinal lift distributions on the cylinders in 
figures 2(a) and 2(b) would be different. 

From the analysis in this section of the wing and cylinder combina- 
tion it can be seen that, for calculating the lift on the cylinder, its 
lateral distribution, and the longitudinal center of pressure, only the 
effect of the vortex potential need be considered. This lift due to the 
vortex potential will be referred to hereinafter as the Induced lift. 

Simplified representation of wing-fuselage combination .- In this 
paper, the wing-fuselage combination will be represented by a system of 
discrete horseshoe vortices and images, so that configurations with 
wings of arbitrary plan forlii may be treated. It is necessary to super- 
pose two pairs of horseshoe vortices of the types shown in figures 1 
and 3 to obtain the vortex- image system of figure k. The vortex- image 
system of figure 4 can be used to represent wings of arbitrary plan 
form in the manner shown in figure 5 l^y locating the bound vortices on 
the wing quarter-chord line or by distributing them over the wing sur- 
face. Since the discrete horseshoe vortex becomes a doublet line eis 
the length of the bo\md leg approaches zero, it is also possible to 
represent the idjig by a continuous distribution of doublets. 

Figure 3 shows the larger pair of vortices with span 2{y-^ + s) 
and strength +F near the smaller pair of span 2(yn - s) and 
strength -F. When the bound legs of the smaller are moved to coincide 
with the bound legs of the larger, the net vortex strength along the 
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section where they coincide is zero, and the remaining sections of the 
larger hound legs form the hound legs of the desired system (see fig. 1+). 
From equations ( 3 ) and (i)-), the lateral distrlhution of the induced lift 
and the total induced lift for the vortex-image system shown in figure If 
eire given hy the following expressions: 



qc 



2(yji - s)^a2 - y2 

(yn - - a2 


tan'l 


in 


2(yn + s)^2 _ y-2 
(yn + - a2 


ALf 

W 



2 a^ 

yn^ - 


( 5 ) 

( 6 ) 


where 



is the loading coefficient on the wing at the station 


yn- 


Since it is known that the longitudinal center of pressure of the 
induced lift for each pair of horseshoe vortices and Images is still 
located at the axis of the hound vortex, it is possible hy superposition 
to calculate the pitching moment on the fuselage in the presence of a 
wing of arhltrary plan form if the lift distrlhution on the wing is 
known, although the conplete longitudinal lift distribution cannot he 
ceilculated unless the conq)onent of loading due to the additional poten- 
tial is calculated. The total Induced lift and its lateral lift dis- 
tribution can also he calculated hy superposition. The method for 
calculating the total induced lift and its lateral distribution and 
longitudinal, center of pressure is discussed in a subsequent section, 
and an Illustrative exan^le is given in appendix B. 


Effect of Finite Fuselage Length on Induced Lift 

In order to obtain an estimate of the error involved in the ass^nI 5 )- 
tlon of the foregoing analysis that the fuselage is infinitely long, an 
approximate calculation will he made of the effect of fineness ratio on 
the Induced lift of an ellipsoid of revolution combined with an infinite 
vortex. The limiting case of the spherical fuseleige is treated first in 
the manner of Vandrey (ref. 5 )^ an*! then the general case is treated. 
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Lift on a combination of sphere and infinite vortex .- The potential 
of the sphere combined -with an infinitf^ vnrt.py (f^g- ) is written 


0 = X tan-1 I + UqX 


(l/x2 + y2 + z2/ 

where the potential for the sphere in the free stream is 


1 + 


(7) 


01 = UqX 


1 + 


Q- 


x2 + y2 + z2 




and the potential for the vortex is 


(8) 


02 = ^tan-l| 


( 9 ) 


The lift on a section of the sphere in the plane y = Constant is 


n\^2.y2 


dLf 


( 10 ) 


where p^ and p^ are the pressiires on the upper and lower surfaces, 
respectively. From Bernoulli's equation, 

P = Po + |p (^ + + (v + Av)^ + (w + Aw )0 (U) 

and 

p^ - p^ = 2 p(u Au + w Aw) 

where Au, Av, and Aw are the velocities induced hy the vortex on the 
upper half of the sphere and u, v, and w are the local velocities on 
the surface of the sphere. Thus 
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Zsa = 


_£ z 

a2 - y2 


Aw = 


_r X 

2 rt a 2 - y 2 


u = 






V = |Uo 


xz 




The lift on the section at y = Constant Is 


dLf 

< 3 y 


3 r pUp 

2 It a^ - y^ 



(12) 


(13) 


which may he written In the form 


since 


dLf 

” P^^ax / ^ 

dy j_f(y) 



(li^) 


and since 


ti — Umax 


when X = 0 , so that 


^ = |Uo 


For an Infinite cylinder with an Infinite vortex^ equation (llf) will also 
hold. In this case 


^nax ~ ^o 
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3,2 _ y2 x2 


and 


f(y) = 00 

From equation (l3) or (lA) the lift obtained for the sphere is 


^ _ 5 

dy 


pUoF 


and for the cylinder 


dLf 

5T - 


Equation (l4) gives the exact resiilt for the case of the sphere 
with an infinite voirtex anfl the same result as the method of reference 5 
for the case of an infinite cylinder with an infinite vortex. There- 
fore, equation (li»-) presumably would give a qualitative estimate for 
the intermediate case of an ellipsoid of revolution having a fineness 
ratio between 1 and infinity. This assumption is made in the following 
section and the value obtained for the induced lift by using equation (l4) 
is shown to be very close to that obtained by the more accurate mathe- 
matical treatment of Vandrey (ref. 5) "tlie case of an ellipsoid 
having a fineness ratio of 

lift on a combination of ellipsoid and infinite vortex .- The local 
lift nn the e lU psoid is given by equation (l^): 

dLf rf(y) 

T — = PUmpv / 2 Au dx 

^ J.fCy) 

From figure 6(b), 



f(y) = d\/l - 


(15) 
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(16) 


The integration indicated in equation (l4) yields 


1 _ ^max/^o 

pUol ^ ^ J 


The values of UmavA^n given as a function of the fineness ratio in 
reference 6 are shown in figure 7* 


1 dLf , 

Figure 8 shows a plot of =- against d/a. The result of a 

pUor dy , 

calculation from reference 5 for the case of an ellipsoid with — = 5 

Q> 

falls near the curve in figure 8. Since the calculations of reference 5 
appear to he accurate to 0 . 05 ^ the agreement may he even better than is 
indicated in figure 8. The results of reference 5 sre obtained by dif- 
ficvilt computations for each body separately, and the method does not 
yield a general expression similar to equation (1T)« 


Figure 8 indicates that the induced lift on the ellipsoids having 

high fineness ratios > 5^ is about 90 percent of the induced lift on 

the infinite cylinder. Since the fuselage is similar in effect to a 

semi-infinite cylinder because of the wake which extends behind it, the 

loss in lift due to the finite fuselage length is about half of that 
Indicated by the calciilation for the ellipsoid of revolution. The veilue 
of induced lift obtained by assuming the fuselage to be an infinite 

^^nax/'^o^ 


cylinder may be multiplied by the factor ^(1 + 


to correct for 


1 + 


the finite fuselage length in the presence of an Infinite wing; however, 
this approximate correction can often be neglected since it is nearly 
unity for most practical fineness ratios. Since the correction is small, 
it is assumed that it may be applied directly for the finite-wing case 

1 L . Ujnax/Uo^ 


without introducing significant error, so that the factor 


+ 


+ ^ 


is to be multiplied by equations ( 5 ) and (6) to correct for finite 
fuselage length. 
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APPLICATION OF SIMPLIFIED MATHEMATICAL MODEL TO CALCULATION 
OF AERODYNAMIC LOADING AND DOWNVJ’ASH 


In the following sections, the method for calculating the induced 
lift will he described and discussed, and methods for calculating the 
spanwise loading on the wing and the downwash will be outlined. The 
methods for calculating the downwash and spanwise lift distribution are 
not rigorous since the effect of the additional, potential, which heis 
not yet been calculated, must be approximated by a simple correction. 
Althou^ the validity of the correction must depend on experimental 
verification, it is believed to be adequate, and the exact value of 
the additional potential may be incoip)orated into the method immediately 
when it is calculated. 


Method for Calculating Induced Lift on Fuselage 

If the lift distribution on the wing in the presence of the fuse- 
lage is known, say from reference 7 (see fig. 9 ) or the method outlined 
in the following section, the induced lift may be calculated very simply. 
For midwing configurations with fuselages of circular cross section, 
equations (5) and (6) are used to calculate the magnitude of the induced 
lift, its lateral distribution, and its longitu din al center of pressure, 
as shown in appendix B. For configurations with arbitrary cross section 
and wing location, the results of appendix A must be used. The method 
of appendix .A, however, will not give the lateral distribution of the 
Induced lift. 

A sanple numerical calculation is shown in appendix B for the mid- 
wing configuration with the circular fuselage. It is necessary only to 
substitute into equations (5) and (6) the value of loading coeffl- 

CCi 

cient for each of the discrete horseshoe vortices . The values of 

c dLf 

the increments of the lateral lift distributions A on the fuselage 

dy 

due to nl 1 of the horseshoe vortices and their images obtained from 
equation (5) are superposed to get the coii5)lete lateral lift distribu- 
tion dLf/dy. The increment of total lift ALf for each pair of horse- 
shoe vortices and images acts at the bound- vortex axis, and the longi- 
tudinal center of pressure of the total lift I^ is obtained simply by 
dividing the sum of the moments of the incremental total lifts AI^ by 
the total lift I^ (the sum of the increments ALf) . 

In calculating the longitudinal center of pressiire of the Induced 
lift, the bound legs of the icing horseshoe vortices should be located 
at the section center of pressure. Since the present methods do not 
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provide a means for calculating the chordwise lift distribution on the 
wing in the presence of the fuselage, the section centers of pressure 
were assumed to lie on the wing quarter-chord line as in approximate 
calculations of the aerodynamic center of the wing alone. However, a 
better approximation to the longitudinal center of pressure of the 
Induced lift may he obtained by assuming that near the wing-fuselage 
juncture and near the wing tip the spanwise variation of the section 
center-of -pressure location is the same as that near the root and tip 
of the wing alone as detennined from measurements or lifting-surface 
calciilations . Curves in reference 8 showing the .spanwise variation of 
section center of pressure for several wing plan forms may serve as a 
guide for estimating the wing-alone values to he iised in preliminaiy 
calculations . 

The results of the calculation for the lateral distribution of 
induced lift and the longitudinal center of pressure canried out in 
appendix B are presented in figures 10 and 11, respectively. Figure 11 
shows the increment of total lift contributed by each p£dr of 

vortices and images and it can be seen that the contribution to the 
total lift Lf due to the outboard part of the wing (the coniponents 
farthest to the right in the figure) is small conipared with the con- 
tribution of the inboard part of the wing. The induced lift on the 
fuselage has also been calculated by using the wing-alone spanwise 
lift distribution shown in figure 9j that is, the effect of the fuse- 
lage on the wing is neglected. It can be seen in figures 10 and 11 
that, although there is about a 10-percent increeiBe in the magnitude 
of dLf/dy and due to the effect of the fuselage on the wing, 

the lateral and longitudinal load distribution on the fuselage is prac- 
tically unaffected. The correction for finite length has not been 
included in the calcvilatlons . This correction would decrease the 
magnitude of the lift, but it would not alter the lateral or longi- 
tudinal distribution. 


The total lift and moment on the fuselage can be obtained by 
adding the components of lift and moment due to the induced lift to 
the components of lift and moment due to the local angle of attack of 
the fuselage (which can be calcvilated as shown in ref. 2). The lift 
and moment on the part of the wing outboard of the fuselage, calculated 
by the method of reference 7 "the method of this paper (described in 
a subsequent section), may then be added to the lift and moment on the 
fuselage to get the total lift and moment on the combination. 
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Outline of Method for Calculating Lift on Wing in 
Presence of Fuselage 

The lift distrihutlon on the wing wiiJL he calculated by equating 
the downwash angle induced at the three-quarter-chord line by the 
horseshoe vortices centered on the quarter-chord line (see fig. 5) to 
the local angle of attank on the three-quarter-chord line at several 
points along the span. Since the boundary conditions on the fuselage 
are not con 5 )letely satisfied by the voitex-image system of figure 5 in 
the region near the bound vortex, it is necessary to resort to certain 
approximations in calculating the downwash. Calculation of the exact 
values of these downwash functions would require the calculation of the 
additional potential, which Involves a great deal of time and effort, 
but if such calculations were made the valties could be used directly 
in the present method and the restidctions suggested in the following 
paragraph would be eliminated. 

The approximations described in the following sections improve as 
the longitudinal distance from the downwash points to the bound vortex 
increases, so that this method is considered to be best suited for con- 
figurations having a fairly small ratio of diameter to root chord, say 

for straight wings with Diafflete r_ ^ 1 swept wings with 

Root chord ^ 

_ Dl _ amete r_ < 1 Calculations of the aerodynamic loading on a wing and 
Root chord 5 

tip-tank conibination, by a method corresponding to the one described 
herein with the downwash points located about 1 tip-tank diameter behind 

the bound vortex j have been found to yield results in 

\Tip chord 

good general agreement with experlmentail results (ref. 9 ). 

Downwash near bound vortex of wing .- Althou^ the boundary condi- 
tions are satisfied conpletely only at infinity, the trailing legs of a 
single peilr of horseshoe vortices and images also satisfy the boundary 
conditions on the cylinder in the plane peipendlcular to the cylinder 
axis which passes through the bound legs of the horseshoe vortices. In 
addition, the boundary conditions on the cylinder are satisfied com- 
pletely by the semi-infinite vortex-image system everywhere in the plane 
of the horseshoe vortex. Since the boxmdary conditions on the cylinder 
are satisfied exactly at the points noted, and are partly satisfied 
everywhere else, it will be assumed that the downwash due to the trailing 
vortices may be calculated approximately, at least in the plane of the 
horseshoe vortex, without introducing any correction factor. 

However, a collection factor must be Tzsed in calculating the do™- 
wash due to the real and image bound vortices which have the greatest 
tendency to violate the boundary condition of zero velocity normal to 
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the surface of the cylinder. In the plane of the wing, along a line 
parallel to the hound vortex, the effect of the cylinder on the vertical 
flow induced hy the hound vortex is assumed to he the same as its effect 
on a two-dimensional uniform rectilinear flow, so that the downwash 
Induced hy the hoimd vortices on that line is increased hy the fac- 

a 2 

tor 1 + 


The downwash angle at a point y = ya x = Xa (shown in 

fig. il-) is then written 


oa = 


1 

WA 2s 



Pn(x,y, 0 ) + Gn(x,y, 0 ) 




(18) 


where the downwash factor due to the trailing vortices and the 
downwash factor due to the hound vortices Gji may he calculated hy 
the Blot-Savart law. 

The function Fq is the sum of four terms, each having the form 




(the term corresponding to the trailing legs of the horseshoe vortex 
centered at y = yh) , and the function Gq Is the sum of four terms, 
each having the form 


1 

1 

yn - ya + 1 

_ _ 


^a " Xn 

/(yn - ya + 1)^ + (xa - Xn)^ 



yn - ya - 1 

- ya - 1 )^ + (xa - 5^)^ 


(the term corresponding to the hound leg of the horseshoe vortex centered 

at y = yn). 
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Scheme for calcvilatinK spanwlse lift distribution on wing .- The lift 
distribution on the -wing is calculated by equating the dovnwash angle on 
the wing induced by the vortex-image system to the local angle of attack 
on the three-quarter-chord line of the wing at m points as in reference 7. 
Thus equation (18) is written for each of the m stations, and the simul- 
taneous equations can then he solved for (-^) m points on the 

span. The + marks in figure 5 Indicate the points where the downwash 
is equated to the local angle of attack. 


The effective angle of attack is equal to the geometric angle 

of attack of the wing Og plus the angle of attack induced by the 
fuselage: 


°-a = °^g + 


and 


1 

Of 


Of 



— (X^ 


= Of 



(19) 


(20) 


where 

tor 


Of is the geometric angle of 


1 + 


takes into account the 


ya 


attack of the fuselage and the fac- 
increase in the vertical velocity 


of the free stream in the neighborhood of the fuselage as calculated by 
assuming the fuselage to be an infinite cylinder in a two-dimensional 
uniform rectilinear flow of magnitude UoOf. 


The effect of finite fuselage length (referred to in ref. 7 as the 

inflow effect ; must be Included separately by multiplying at 

c 

each spanwise station by the f ant or 1 + 2e, where e is the ratio of 
the local Increment of longitudinal velocity due to the fuselage to the 
free-stream velocity. As shown in reference J, the factor 1 + 2 e is 
\ised to account approximately for the small Increase in dynamic pressiire 
of the flow over the wing due to the increase in the local longitudinal 
velocity near the surface of a fuselage of finite length. 
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Outline of Method for Calculating Downweish 

The downwash calculations made in this section are for the case of 
the wing at an angle of incidence with the fuselage at zero geometric 
angle of attack. The accuracy of the calculated results must he verified 
by experimental results; however, by coEoparison with results of doirawash 
calculations for the wing alone, the results of the calculations for the 
downwash behind the wing-fuselage combination may serve to give useful 
information regarding the effect of the wing-fuselage Interference. 

As was noted previously, the boundary conditions on the wing-fuselage 
combination are conpletely satisfied by the vortex- image system infinitely 
far behind the wing. At a great distance behind the wing, therefore, 
the theoretical value of the downwash can be calculated accurately for 
a given spanwise lift distribution. This calculation is simple, and the 
resxilts may be useful for certain applications. However, the approxima- 
tions made in the preceding section for calculating the downwash at the 
wing three-quarter-chord line may give more accurate results for the 
downwash in the plane of the wing in the region nearer the bound vortex. 

If the lift distribution on the wing is known, say from the method 
of reference 7 or "tlie method outlined in this paper, the downwash 
in the plane x = 00 may be calculated by adding the downwash due to 
each pair of trailing vortices and their images. Thus, from the Biot- 
Savart law and figure k, the downwash angle Oa at the point y = ya, 
z = Za, and x = “ is given els 

(.1) 

where 


Fn = 


^n - ya + 1 


yn - ya - 1 


yn + ya + 1 


(yn - ya + + 2^2 (y^ - ya - l)^ + (yn + ya + l)^ + Za^ 


yn + ya - 1 


ya - ^ + 1 
7n 


- -S2 

ya - — - 1 
yn 


(yn + ya - 1)^ + ^ 


a^ 


- 2 


ya - ^ - 1 + za‘ 

^ yn / 


— a^ 


_ a2 




(22) 


and (cc;/c)j^ is the loading coefficient at station iy^. 
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The representation of the wake hy discrete vortices is satisfactory 
if the downwash is calculated at points halfway between the trailing 
legs (yg, = yn) • The values calculated at these points may he faired to 
obtain a continuous spanwise distribution of downweish angle. 


COWCLUDING REMARKS 


For the purpose of calciilating the longitudinal loading on the 
fuselage in subsonic flow, a midwing wing-fuselage combination with a 
fuselage of circular cross section has been represented by a system of 
discrete horseshoe vortices and images. By using this simplified 
mathematical model, or the extension of it given in an appendix for 
nonmldwing configurations with ftiselages of arbitrary cross section, a 
method is derived for cadcxilating the lift on the fuselage induced by 
the wing. The method is i3JLustrated by a numerical exaniple. This 
"induced lift' can be added to the 31ft on the part of the wing out- 
board of the fuselage and the 3dft on the fuselage due to the upwash 
induced by the wing to get the total loading on the wing-fuselage com- 
bination. 

In addition to the method for calculating the induced lift, which 
is theoretically rigorous, methods for calculating the downwash far 
behind the wing and for calculating the spanwise lift distribution on 
the wing for midwing configurations with axisymmetrlc fuselages ere 
outlined. In calculating the spanwise lift distribution on the wing, 
approximations are made to account for the effect of the "additional 
potential, so that the method is not rigorous. However, the effect of 
the additional potential may easily be Incorporated into the method 
when it is calculated. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 15, 1952- 
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APPEiroiK A 

CALCULATION OF INDUCED LIFT FOR CONFIGURATIONS WITH 
ARBITRAEOT CROSS SECTION AND WING LOCATION 


This appendix indicates how the relation, given in reference 1)-, 
between the magnitude of the Induced lift on a cylindrical fuselage 
and the lift distrihutlon on the wing can he generalized so that a 
siD 5 )le equation can he obtained from which the longitudinal center of 
pressaire of the Induced lift may he calculated. Reference 4 shows that 
the following relation will hold between the total lift induced on a 
cylindrical fuselage of arbitrary cross section an(i the lift distribu- 
tion on a wing in the horizontal plane z = Constant: 


= / pUor(y)g(y)dy (A1) 

Jw 

where y is the lateral ordinate, I*(y) is the spanwise distribution 
of circulation, g(y) is the increment in the flow velocity normal to 
the plane of the wing produced by a unit uniform transverse flow normal 
to the cylinder axis, and W indicates that the Integral is taken over 
the span of the exposed wing. 

In order to generalize equation (A1) so that the center of pressure 
of the induced lift on the fuselage can be calculated, it is convenient 
to consider wing elements to be represented by horseshoe vortices having 
bound legs of vanishingly small length 5y^ lying in the horizontal 

plane z = Constant with trailing legs parallel to the x-axls (longi- 
tu d i n al axis). On a cylinder having its axis coinciding with the x-axls, 
these horseshoe vortices will each induce a lift of magnitude 5Lp, 
defined by 


^n+5yn 

Sir = J pUor(y)g(y)dy (A2) 

yn-syn 

where is the lateral location of the center of the bound leg. It 

can be shown that the longitijdlnal location of the center of pressure 
of this induced lift 61^ is at x^, the same longitudinal location as 
the bound leg of the horseshoe vortex. The proof is exactly the same in 
all essentials as that presented in the body of this paper for the case 
of a circular cylinder with a finite horseshoe vortex located in the 
plane z = 0. 
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From ttiese considerations, it is clear ttiat the total induced lift 
on the fuselage and its longitudinal center of pressure can he 

given as follows; 


Lf = 




pUo / 7 (x,y)g(y)dS 
-'S 

pUo X 7 (x,y)g(y)dS 


(A5a) 


(A5b) 


where 7 (x,y) is the strength of an infinitesimal horseshoe vortex 
centered at (x,y), and the integrals are to he taken over the exposed 
■\7ing STirface S. It may he noted that, -vdiereas equation (A5a) is valid 
at all suhsonlc and supersonic Mach nunibers, equation (A^b) is valid only 
at subsonic Mach numbers- 

Equations (A5a) and (AJb) represent, in general, the desired equa- 
tions for the magnitiide and center of pressure of the induced lift on a 
cylinder in the presence of a lifting surface. In order to illustrate 
their application, they will he developed for the case of a circular 
fuselage with a lifting line of horseshoe vortices and will he reduced 
so that the siitple computational procedure of appendix B can he used. 

For a circular cylinder of radi\is a, g(y) is given as 

g(y) = =^ 1 ) . 

(y2 + z2)2 

When the wing is represented by a lifting line, equation (A5a) has the 
same form as equation (A2) and the increment in lift ALf due to two 
horseshoe vortices of strength r having spans 2s and located in the 
horizontal plane z = zji at +yn and -y^^ is 

Alf = 2pUor a2(y2 - zf,2) ^ 

yn-s (y2 + Zj^2)^ 


= a2s* 

zji^ + 2zii2(yj^2 + b2) + (y^2 - 


W 


It may he noted that for the midilng case zj^ = 0, equation (a 4) is 
identical to equation (6). 
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APPENDIX B 

imJSTRATIVE EXAMPIE OF INDUCED- LIFT CALCULATION 


A numerical exan^jle is given to illustrate the method presented in 
the hody of the paper for calciilating the Induced lift on the fuselage. 


Geometric Characteristics of Configuration Used 
in Illustrative Exanple 

The plan view of the wing-fuselage comhinatlon is essentially the 
same as shown in figure 5 « geometric data are: 


Aspect ratio 8 

Taper ratio 0.^5 

Sweephack, deg 45 

a* 0.10 

s* 0.05 


Spanwlse Loading on Wing in Presence of Fuselage 

The spamd-se lift distrihution on the wing in the presence of the 
fuselage for a hl^-midwing configuration, tahxilated below and shown in 
figure 9 , was obtained by the method of reference 7. It is qualitatively 
correct for the pure midwing case and is used to illustrate the procedure 
for obtaining the Induced lift on the fuselage and its lateral distribu- 
tion and longitudinaJ. center of pressure. The lift distribution on the 
wing in the presence of the fuselage, plotted in figure 9 > Is tabulated 
as follows : 


n 

yn* 

xn* 

(-) 

\ c /n 

1 

0.15 

0.15 

0.569 

2 

•25 

•25 

.566 

3 

•55 

•35 

• 556 

4 

•45 

•45 

.558 

5 

•55 

•55 

• 519 

6 

.65 

• 65 

.300 

7 

•75 

•75 

.266 

8 

.85 

.85 

.241 

9 

•95 

•95 

.200 




2h 


MCA TN 5057 


Longitudinal Distribution of Induced Lift 

Equation ( 6 ), Tdiicb gives the lift on the fuselage due to a single 
pair of horseshoe vortices at (xn',±yn') and their ijnages, can be 
written 


^ S. ^ 

<13 \ c /n ~ ® 

The total lift 1 ^/q.S is then 


q.S 



(Bl) 


(B 2 ) 


The longitudinal center of pressure referred to the intersection of the 
quarter-chord line and the fuselage axis of symmetry (see fig. 5 ) is 


Xcp 


The coD5)uted values of Lf/q.S and x^p* are obtained from the values 
presented in the following table; 



n 

2 

RBjjHjHH 


1 2 1 2 

yn - s 


U/nyn’ 2 - b ’2 

1 

1.000 

0.369 

0.55^ 

2 

.533 

.122 

.305 

3 

.168 

.060 

.210 

k 

.100 

.054 

-.153 

5 

.066 

.021 

.116 

6 

.dks 

. 0 ll<- 

.091 

T 

.036 

.010 

.075 

8 

.028 

.007 

.060 

9 

.022 

.005 

.048 

■ 


= 0 . 614-2 

= 1.612' 
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qS 


0.05 


fcc-. 


C /ny^*2 . g.2 


0.032 


X 


cp 


I 



2.51 


The values of 


Wny„’2- s'2 


are plotted as vectors in 


figure 11. The location of x^,p' is also shown. 


Lateral Distrihution of Induced Lift 

Equation ( 5 ), which gives the lateral distribution of the induced 
lift on the fuselage due to a single pair of horseshoe vortices at 
(xn'^iyn') and their images, can be written 




2(yn' - s' )jl - y'^ 
(yn' - s’)^ - 1 


tan"^ 


2 ty„' + b')/i - y'g 
+ B'f - 1 


(Bit) 


The total lift at a station ±y' on the fuselage is then 




taE-l 

1 

m 

10 

. tan-l 

+ e')\ll - y'^ 

qc ^ 

V a / n 

(yn' 

- B'f - 1 

(yn' 

' + s’)^ - 1 


(B5) 
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The value of 


<U^/dy 


at ±y’ = 0.25 


q.c 

the following table (the lateral lift 
lated in a similar manner); 


is 

at 


calculated from the values in 
any station ±y' can he calcu- 


1 

© 




_ 2(y^' - b')/i - y’2 

2(yn' + s')\/l - y’2 

((^ - (2)') 

n 

tan o 

(yn - b')^ - 1 

(y„- s')2 - 1 

V c 

1 

1.571 

0.918 

0.241 

2 

.918 

.653 

.104 

5 

.653 

.481 

.054 


.481 

.387 

.052 

5 

.387 

.324 

.020 

6 

.524 

.279 

.014 

7 

.279 

.244 

.010 

8 

.244 

.217 

.007 

9 

.217 

.195 

.005 

1 



= 0.487 


Thus, at y' = 0.25, 


/ 


= |(0.1^8T) = 0.510 
qc 


The con5)lete lateral lift distribution is shown in figure 10. 
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(b) Infinite-vortex — ellipsoid configuration. 


Figure 6.- Illustrations of sphere with infinite vortex and ellipsoid 

with infinite vortex. 




~ - Asymptote at Infinity — 
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lution. 
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Rctio of local lift on ellipsoid to local Ii1^ on 
. 'dLf 



Figure 8.- Effect of fineness ratio on the lift on an ellipaoid of 

revolution intersected by an Infinite vortex having the strength P . 
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Calculated with fuselage effect on wing 
Calculated without fuselage effect on wing 



.2 .3 ,4 .5 .6 .7 .8 .9 UO 


Ointensioniess lateral ordinate, y' 

Figure 10.- Letereil distrlljutlon of indoiced lift on the fuselage in the 
presence of the vlng. Sweephack angle, 45°; taper ratio, 0.45; 

8} a* => 0.10. 
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Dimensionless Increment of lift 


CD 



/-Wing center section 
— Root section 


Tip section 




4- 


I 23456789 
Dimensionless longitudind ordinate, x’ 


10 


Figure 11.- Longitudinal distribution of induced lift on the fuselage in 
the presence of the wing. Sweephack angle, taper ratio, 0.55? 

A “ 8; a* =. 0.10. 
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